Focusing of light by random scattering 
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Random scattering of light is what makes materials such as white paint, clouds and biological 
tissue opaque. We show that although light propagating in these media is diffuse, a high degree of 
control is possible as phase information is not irreversibly lost. Opaque objects such as eggshell or 
white paint focus coherent light as sharply as a lens when illuminated with a wavefront that inverts 
the wave diffusion. We demonstrate the construction of such wavefronts using feedback, achieving 
a focus that is 1000 times brighter than the diffusely transmitted light. Our results are explained 
quantitatively by a universal relation based on statistical optics. 

PACS numbers: 87.64.Cc, 42.30.Ms, 42.25.Dd 
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Optical microscopy and spectroscopy are essential 
tools for the study of living organisms and inanimate 
objects. These methods rely on the ability to deliver and 
collect light with a high degree of control. Unfortunately, 
in many organic and inorganic materials light is scattered 
randomly and the required directionality of the light is 
lost P, 0, d- I n these materials, light performs a ran- 
dom walk and emerges as a random interference pattern 
known as speckle. Innovative imaging methods are di- 
rected towards isolating the unscattered fraction of the 
light jj,0| and towards obtaining useful information from 
the multiply scattered light 0, 0, H • Ideally, one would 
completely eliminate or counteract scattering. 

We demonstrate inverse wave diffusion, a method for 
counteracting scattering and diffusion of light. By con- 
structing a perfectly matched wavefront, we make nor- 
mally opaque objects focus light as sharply as a lens. This 
wavefront cannot be known a priori and is constructed 
using feedback from a detector in the target focus. By 
changing the incoming wavefront, we control the position 
and shape of the focus; it is even possible to transmit col- 
limated beams or simple images. Inverse wave diffusion is 
universally applicable to scattering objects regardless of 
their constitution and scattering strength. We envision 
that, with such active control, random scattering will be- 
come beneficial, rather than detrimental, to i mag ing [9j, 
communication [lfjl Hlj and non- linear optics |l3 | . 

Figure ^ shows the principle of the experiment. Nor- 
mally incident light from a 632.8 nm HeNe laser is scat- 
tered by the sample and forms a random speckle pattern 
(Fig. [fla)). The goal is to match the incident wavefront 
to the sample, so that the scattered light is focused in a 
specified target area (Fig. njb)). We divide the incident 
wavefront into an array of spatial segments. Initially, 
the amplitude from each segment contributes randomly 
to the light amplitude in the target focus (Fig. [He)). 
With a spatial phase modulator, the wavefront of the in- 
cident light is adjusted using the intensity in the focus 
as feedback. The phase of each segment is adjusted until 
the target intensity is maximal. After adjustment, the 
phase modulator inverts the diffusion in the sample and 
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FIG. 1: Design of the experiment, (a) A plane wave is focused 
on a disordered medium, a speckle pattern is transmitted, (b) 
The wavefront of the incident light is shaped so that scatter- 
ing makes the light focus at any desired point, (c) Complex 
amplitude representation of the field at the target. Before op- 
timization, each segment of the incident wavefront contributes 
to the field (thick arrow) in a random way. By adjusting the 
phase of a single segment, we determine the phase at which 
the total field is maximal (dotted arrows), (d) The phase of 
the incident light is adjusted to have all fields interfere con- 
structively. The target intensity is at the global maximum. 



the scattered light interferes constructively in the target 
focus (Fig. Old)). 

We performed first tests of inverse wave diffusion using 
rutile TiC>2 pigment, which is one of the most strongly 
scattering materials known. The sample consists of an 
opaque, lO.T/mi thick layer of white pigment |13| with 
a transport mean free path of 0.55 ± 0.10/im measured 
at a wavelength of 632.8 nm (See Appendix A). Since in 
this sample the transmitted light is scattered hundreds 
of times, there is no direct relation between the incident 
wavefront and the transmitted image [lj, [I3 • Figure [21 
shows the intensity of the transmitted light seen through 
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FIG. 2: Shaped transmission through a strongly scattering 
sample consisting of Ti02 pigment, (a) Transmission micro- 
graph with an unshaped incident beam. The scattered light 
forms a random spackle pattern, (b) Transmission after opti- 
mization for focusing at a single target. The scattered light is 
focused to a diffraction limited spot that is 1000 times brighter 
than the original speckle pattern (c) Multi-beam optimiza- 
tion. The disordered medium generates five sharp foci at the 
defined positions. FiguresEJa) toEIc) are represented on the 
same logarithmic color scale that is normalized to the average 
transmission before optimization, (d) Correlation function of 
the incident wavefront used to form the pentagon seen in Fig. 
|2fc). There is no correlation between neighboring segments, 
which indicates that their contributions to the target function 
are independent. The incident wavefronts are composed of a 
total of 3228 individually controlled segments. 



a microscope objective. In the first image (Fig. Eta)) we 
see the pattern that was recorded when a plane wave was 
focused onto the sample. The transmitted light formed a 
typical random speckle pattern with a low intensity. We 
then optimized the wavefront so that the scattered light 
focused to a target area with the size of a single speckle. 
The result is seen in Fig. ED 3 ), where a single bright spot 
stands out clearly against the diffuse background. The 
focus was over a factor 1000 more intense than the non- 
optimized speckle pattern. Instead of focusing light to a 
single spot, it is also possible to optimize multiple foci 
simultaneously. By adjusting the target function used 
as feedback, the scattered light was made to form sim- 
ple images, such as the pentagon shown in Fig. Et c )- 
Figure Eld) shows that there is no correlation between 
neighboring segments of the optimal incident wavefront, 
which indicates that the light was fully scattered and no 
ballistic transmission occurred. 

The TiC>2 produced a high quality focus, as can be 
seen in Figure EJ For comparison, the intensity auto- 
correlation function of the speckle before optimization is 
also shown. From statistical optics it is known that the 
profile of the autocorrelation function equals the diffrac- 
tion limited beam profile |15f. Since the two functions 
overlap, we conclude that the multiply scattering TiC>2 
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FIG. 3: Cut through the profile of the target focus after op- 
timization (see Fig. Et^)) compared to the normalized auto- 
correlation function of the speckle before optimization (dotted 
line) . 



Sample 


L (/jm) 


d (/im) 


f (mm) 


r, (±15%) 


N 


Ti0 2 


10.1±0.3 


18.7 


3.5±0.5 


1080 


3228 


Petal, fresh 


43±5 


10.6 


3.5±0.5 


64 


208 


Petal, dry 


37±5 


12.6 


3.5±0.5 


630 


1664 


Eggshell 


430±30 


2.1 


3.5±0.5 


250 


3228 


Tooth 


1500±100 


155 


125±10 


70 


208 



TABLE I: Intensity enhancement for different materials. L, 
sample thickness (± surface roughness), d, diameter (1/e 2 in- 
tensity) of focus, /, distance between random medium and 
focus, r\, maximum enhancement reached, N, number of seg- 
ments used by the algorithm to describe the wavefront. 



produces a diffraction limited focus, i.e. the disordered 
layer of pigment focuses light as sharply as an ideal lens 
of the same size. These results are clear proof that the 
wavefront of multiply scattered transmitted light can be 
controlled with high accuracy. 

Inversion of wave diffusion was realized in a variety of 
opaque materials ranging from white pigment to fresh 
flower petals. Table [I] summarizes the results for the 
different materials used. Although the samples vary 
in thickness, composition and scattering strength, they 
were all able to focus a properly prepared wavefront to a 
diffraction limited spot. The intensity enhancement r\ - 
defined as the ratio between the optimized intensity and 
the average intensity before optimization - varies between 
60 and 1000. The main reason for this variation is that 
the temporal stability of the transmitted speckle pattern 
is not the same for all materials. 

The optimization procedure makes use of the linear- 
ity of the scattering process. The transmitted field in 
the target area, E m , is a linear combination of the fields 



coming from the N different segments of the modulator, 
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71 = 1 
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(1) 



where A n and <j> n are, respectively, the amplitude and 
phase of the light reflected from segment n. All scattering 
in the sample is described by the elements t mn of the 
unknown transmission matrix. Clearly, the magnitude 
of E m will be the highest when all terms in Eq. ^ are 
in phase (also see Fig. ^c), IHd)). We determine the 
optimal phase for a single segment at a time by cycling its 
phase from to 2tt. When the phase of a single segment 
n is changed, the target intensity detected by the CCD 
responds as 
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\E m \ — \Eq\ + \t mn A n \ + 



2\E \\t mn A n \ cos [axg(t mn ) - aig(E ) 



(2) 



where Eq is the field of the scattered light originating 
from all segments except segment n. When N is large, 
each segment contributes little to the total field and Eq 
is equal for all segments. For each segment we store the 
phase at which the target intensity is the highest. At 
that point the contribution of segment n is in phase with 
the already present diffuse background Eq. After the 
measurements have been performed for all segments, the 
phase of the segments is set to their stored values. Now 
the contributions from all segments interfere construc- 
tively and the target intensity is at the global maximum. 
This method is generally applicable to linear systems and 
does not rely on time reversal symmetry or absence of ab- 
sorption. 

The maximum intensity enhancement that can be 
reached is related to the number of segments that are 
used to describe the incident wavefront. For a disor- 
dered medium the constants t mn are statistically in- 
dependent and obey a circular Gaussian distribution 
|15l Ha. ll7L HH ligj and the expected enhancement 77 can 
be calculated, 
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It was assumed that all segments of the phase modula- 
tor contribute equally to the total incident intensity. We 
expect the linear scaling behavior to be universal as Eq. 
121 contains no parameters. Neither sample thickness nor 
scattering parameters will influence the expected inten- 
sity enhancement. Also, since we are free to choose the 
basis for Eq. Q] we expect to find the same enhancement 
regardless of whether the target is a focus or a far-field 
beam and regardless of how the shaped wavefront is pro- 
jected onto the sample; of course the required optimal 
wavefront will be different for these varying configura- 
tions. The number of degrees of freedom N is bound to 



FIG. 4: Enhancement as a function of the number of seg- 
ments. The experiment was performed twice, once with the 
sample in focus (squares) and once with the sample 100 /im 
out of focus (triangles). The solid line represents the theo- 
retical enhancement for an ideal system (Eq. 0. The dot- 
ted line represents the expected enhancement when residual 
amplitude modulation of the phase modulator and the finite 
persistence time of the sample are taken into account. The 
uncertainty in 77 is of the order of the symbol size. 



a maximum given by the number independent speckles 
on the sample surface; N max = 8A/X 2 , where A is the 
illuminated surface area of the sample and A is the wave- 
length of the light |13)lly|- When TV approaches N max , 
the intensity in the target focus alone becomes compara- 
ble to the total transmission before optimization. In this 
extremely interesting regime, the assumptions underlying 
Eq. |21 arc no longer valid. With our current apparatus 
N <C iV max . 

We tested the universal scaling behavior implied by 
Eq. |21 by changing the number of segments into which 
the phase modulator is subdivided. Using the same Ti02 
sample as before, the algorithm was targeted to construct 
a collimated beam. In Fig. 01 the enhancement is plot- 
ted as a function of the number of segments for different 
focusing conditions. The linear relation between the en- 
hancement and the number of segments is evident until 
the enhancement saturates at 77 = 1000. All measured 
enhancements were slightly below the theoretical maxi- 
mum. This is understandable since the phase space is 
huge and all perturbations move the system away from 
the global maximum. The main reason for deviations 
from the optimal wavefront is residual amplitude mod- 
ulation in the phase modulator. The amplitude modu- 
lation introduced a constant, uncontrollable bias in the 
field (See Appendix A). 

The saturation of the enhancement was the result of 
slow changes in the speckle pattern. This instability ef- 
fectively limited the number of segments for which the 
optimal phase could be measured. We estimate that the 



effective enhancement decreases to rj s = r]/(l + NT/T s ), 
where T = 1.2s is the time needed for one measurement 
and T s = 5400s is the timescale at which the speckle pat- 
tern of the TiC>2 sample remains stable. The instability 
of the speckle pattern was most likely caused by a fluc- 
tuating humidity of the sample [20]. Depending on the 
environmental conditions, T s can be considerably higher 
and enhancements of over two thousand have been mea- 
sured overnight. 

Our results show that precise control of diffuse light is 
possible using an optimal, non-iterative algorithm; light 
can be directed through opaque objects to form one or 
multiple foci. The brightness of the focal spot is ex- 
plained by a model based on statistical optics. We expect 
inverse wave diffusion to have applications in imaging 
and light delivery in scattering media. Dynamic mea- 
surements in biological tissue are possible when the time 
required for achieving a focus can be reduced to below 1 
ms per segment jj, |(j ; we estimate that this timescale is 
technologically possible with the use of fast phase modu- 
lators. Furthermore, the high degree of control over the 
scattered light should permit experimental verification 
of random matrix theories for the transport of light and 
quantum particles [H, [Hi ■ 
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tensity coming from the phase modulator is 1.0 2 /(0.4 2 + 
1.0 2 ) = 0.86. After each adjustment of the phase, we al- 
low for a stabilization time of 100 ms. The phase shaped 
beam is spatially filtered to remove higher order diffrac- 
tion and focused onto the sample using an objective with 
a magnification of 63 times and a numerical aperture of 
0.85. Since the solid angle covered by the objective is 
0.957T and only a single polarization is used, at most 24% 
of the mesoscopic channels can be addressed. A fraction 
of the phase shaped light is directed to a silicon photo- 
diode which acts as an intensity reference. 



Detection 

The transmitted light is imaged using a microscope ob- 
jective with an NA of 0.5 and a magnification of 20 times. 
The light passes a Glan Thompson polarizer and is de- 
tected using a 2/3" , 12-bit CCD camera in the back focal 
plane of the objective. The camera image is integrated 
over an area that is smaller than the typical speckle size. 
No collimating optics are used behind the sample. Per 
segment of the phase modulator, the phase is varied be- 
tween and 2n in ten equal steps and a sinusoid is fitted 
to the measured intensities. To minimize the effect of 
noise, the optimization procedure first performs a coarse 
pre-optimization using 12 segments. Then, the algorithm 
is executed twice. The result of the first iteration is used 
as the reference field for the second iteration. The av- 
erage intensity of the speckle background is obtained by 
averaging over 4000 random phase patterns. 



Samples 



I. APPENDIX A: MATERIALS AND METHODS 

Phase shaping 

The experiments are performed using a polarized 5 
mW Helium-Neon laser with a wavelength of 632.8 nm. 
Phase modulation is achieved using a Holoeye LCR-2500 
twisted nematic liquid crystal reflective spatial light mod- 
ulator. The experimental configuration is shown in Fig- 
ure Al. We illuminate a circular area containing 3 ■ 10 5 
pixels grouped together in square segments. The phase 
modulator operates in a phase mostly mode[2l|. In this 
mode, the field modulation curve can be closely approx- 
imated by a circle in the complex plane. Due to residual 
amplitude modulation, the centre of this circle is biased 
with respect to the origin. The bias corresponds to an 
uncontrollable fraction of 40% of the field. Therefore, the 
ratio of the controlled intensity to the average total in- 



Four different samples were used: Ti02, flower petals, 
egg shell, and a primary tooth. The first sample consists 
of an opaque, 10.1 /zm thick layer of rutile Ti02 pig- 
ment on a 2mm-thick fused silica substrate lj- By mea- 
suring the total transmission, the transport mean free 
path was found to be 0.55 ± 0.10/im at a wavelength 
of 632.8 nm. The sample was placed with the pigment 
layer towards the first microscope objective. The flower 
petal was freshly picked from a Bellis Perennis (Common 
Daisy) and fixated with room temperature parafilm be- 
tween a microscope slide and a coverslip. A second petal 
was wet mounted on the slide after which the sample was 
allowed to dry for one day. The egg shell is from a white 
chickens egg. A part of the shell was rinsed, dried and 
placed between the microscope objectives. The primary 
tooth (incisor) was placed in the focal plane of the first 
microscope objective. In this latter case, a second objec- 
tive was not used. 
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Figure Al: Schematic drawing of the apparatus. 
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